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Abstract. We study the intermediate statistics of the spectrum of quasi- 
energies and of the eigenfunctions in the kicked rotator, in the case when 
the corresponding system is fully chaotic while quantally localized. As 
for the eigenphases, we find clear evidence that the spectral statistics 
is well described by the Brody distribution, notably better than by the 
Izrailev's one, which has been proposed and used broadly to describe 
such cases. We also studied the eigenfunctions of the Floquet operator 
and their localization. We show the existence of a scaling law between the 
repulsion parameter with relative localization length, but only as a first 
order approximation, since another parameter plays a role. We believe 
and have evidence that a similar analysis applies in time-independent 
Hamilton systems. 



1 Introduction 

One of the most important manifestations of quantum chaos of low-dimensional 
classically fully chaotic (ergodic) Hamiltonian systems is the fact that in the 
(sufficiently deep) semiclassical limit the statistical properties of the discrete 
energy spectra obey the statistics of Gaussian Random Matrix Theory (RMT). 
The opposite extreme are classically integrable systems, which quantally exhibit 
Poissonian spectral statistics (see pQ). 

Quantum kicked rotator (QKR) is a typical example in the field of quantum 
chaos 2\ . A typical property of the QKR is the chaos suppression for sufficiently 
large time scales. The study of the statistical properties of the classical and 
quantum (semiclassical) parameters in such systems is of great importance. Here 
we study in detail the semiclassical region where k > K > 1, i.e. the regime of 
full correspondence between quantum and classical diffusion (on the finite time 
scale t < t]j) and the manifested quantum dynamical localization for t > to- 
Furthermore, we are focused in the probability level spacing distributions in the 
regime where the system is classically strongly chaotic (K > 7) but quantally 
localized, i.e in the intermediate or soft quantum chaos, as it is described in the 
literature [51. 
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2 The quantum kicked rotator model 



The QKR model [3] is described by the following function 
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where ft is Planck's constant, / is the moment of inertia of the pendulum and 
£o is the perturbation strength. The motion after one period T of the ip wave 
function then can be described by the following mapping 



ip(6,t + T) = Uip{0,t), 
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where the ip function is determined in the middle of the rotation, between two 
successive kicks. The evolution operator U of the system corresponds to one 
period. Due to the instant action of the perturbation, this evolution can be 
written as the product of three non-commuting unitary operators, the first and 
third of which corresponds to the free rotation during half a period G(r/2) = 

exp^i^j^j^, t = tiT/I, while the second B(k) = exp(— ik cos 9), k = e /h 
describes the kick. The system's behavior depends only on two parameters, i.e. 
r and k and its correspondence with the classical systems is described by the 
relation K = kr = e T/I. In the case K = kr >> 1 the motion is well-known 
to be strongly chaotic. The transition to classical mechanics is described by the 
limit k — > oo, r — > while K = const. In what follows h = t and T = I = 1. We 
shall consider mostly the semiclassical regime k > K, where r < 1. 

In order to study how the localization affects the statistical properties of the 
quasienergy spectra we use the model's representation with a finite number N 
of levels [315] 
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m— 1 



The finite unitary matrix U nm determines the evolution of a N- dimensional 
vector (Fourier transform of ip) of the model 
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where Gu> — exp [irl 2 /4J Sw is a diagonal matrix corresponding to free rotation 
during a half period T/2 and the matrix B n > m i describing the one kick has the 
following form 
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The model (j4]) with a finite number of states is considered as the quantum 
analogue of the classical standard mapping on the torus with closed momentum 
p and phase 8 where U mn describes only the odd states of the systems, i.e. 

3 Intermediate statistics and comparison of probability 
distributions 

Let us first compare the Brody and Izrailev probability distribution functions 
(PDFs) for the study of the intermediate level statistics. The Brody distribution 
is defined by the relation 



where the two parameters C\ and C% are determined by the normalization con- 
ditions J °° Pbr{s)(Is = 1 and f °° sPb(s)cIs = 1. Izrailev suggested the following 
distribution (see [5] and references there for the details and the argumentation) 



in order to describe the intermediate statistics, where the parameters A and B 
are determined again by the two above normalization conditions. Both /3 param- 
eters, in the strongly localized regime tend asymptotically to with Poissonian 
statistics while in the chaotic one tend to 1, which excellently approximates 
the Gaussian Orthogonal Ensemble (GEO) of the RMT. On the other hand, 
the non-integer /3 in the PDFs could be associated with the statistics of the 
quasienergy states with chaotic localized eigenfunctions. Here, we use N — 4000 
(which is considerably much larger size compared to the one used in the past 
studies) and K = 7 with k — 30. In Fig. [IJa) we show the numerical data (his- 
togram) and the two PDFs. Their repulsion parameters have been calculated 
with best fit procedure independently. The corresponding values are found to be 
/?br ~ 0.424 and (3\z ~ 0.419 for the two PDFs respectively. In the inset figure 
of Fig.QJa), we may see how the -Pbr(s) manages to capture and describe better 
the peak of the distribution where the most significant part of quasienergies uj 
is concentrated. The dot-dashed gray line indicates the Wigner surmise while 
the dot-dot-dashed one the Poisson distribution. Similar findings have also been 
found even for smaller sizes of the matrix U nm , where the statistics are improved 
by sampling more matrices with slightly different values of k as e.g. in [5]. 

The above statement, regarding the Pbr(s) better agreement with the numer- 
ical data, becomes more clear when checking the so-called [/-functions U(s) = 
(2/tt) arccos ^/l — W(s) of the two above distributions [6]. The W(s) = P(x)dx 
is the cumulative (or integrated) level spacing distribution function (CDF). The 
[/-function has the advantage that its expected statistical error SU is indepen- 
dent of s, being constant for each s and equal to SU — l/(n\/N s ), where N s is 
the total number of objects in the W(s) distribution. The numerical pre-factor 
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Fig. 1. Intermediate statistics (panel (a)) for distribution P(s) (histogram - black solid 
line) of the model (H}-([5} fitted with distribution Pbr(s) (black dashed line) and Piz(s) 
(black dotted line) for N = 4000, K = 7 and k = 30 (see text for discussion). The gray 
lines indicate the two extreme distributions, i.e. the Poisson and Wigner. In panels 
(b),(c),(d) we show the comparison of the the Brody (black line) and Izrailev (gray 
line) PDFs with the numerical using the [/-function and ^-distribution (see text for 
discussion). 



2/ir is determined in such a way that U(s) e [0, 1] when W(s) € [0, 1]. We may 
note here that the (3 values for the CDFs may be in principle slightly different 
compared to those found by the PDFs. In Fig. QJb), we show the Ubr — U and 
Uiz — U vs. Ubr where one may see that the Brody one is in general closer 
to zero (black line) than the Izrailev one (gray color). This fact indicates that 
the Brody one fits better the numerical data. This is also evident in Fig. [He), 
where the Wbr — W and W\z — W vs. Wbr are presented while in Fig. QJd) 
the Wbr, — W and W\z — W vs. s. The horizontal zero line in these panels 
indicates the complete agreement between the numerical data and theoretical 
predictions. The repulsion parameters for the CDFs used in panels (b),(c),(d) 
are /3br ~ 0.396 and /3iz ~ 0.366 respectively. 

4 Scaling laws and localized chaotic regimes 

A number of different ways to measure and estimate the localization length of 
the eigenfunctions have been proposed in the literature. Here, we adopt the well- 
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accepted measure described and justified in e.g. [5]: For each iV-dimensional 
eigenvector of the matrix U nm the information entropy is 3%n(ux, Ujsr) = 
~ 5Zn=i u n ^ n u ni where u n = Re ip n and J2n u n = 1> The distribution of it„ 
for the GOE in the large iV-limit tends to the Gaussian distribution and we get 
jjfGOE = ^ (0.5JV + 1) -ip(l.5) ~ ln(0.5iVa) + O(l/iV), where a = 4/exp(2- 
7) ~ 0.96, while ^ is the digamma function and 7 the Euler constant (~ 0.57). 
Then the entropy localization length Ih is defined as Ih = N exp (J#n — Jf^ OE ^ . 
The fluctuations can be minimized when using the mean localization length < 
Ih >= d, which is computed by averaging over all eigenvectors of the same 
matrix (or over an ensemble of similar matrices) d = N exp (< J^m > —J^n° E )- 
The parameter that determines the transition from weak to strong quantum 
chaos is not the strength parameter k but the ratio of the localization length Zoo 
to the size N of the system, A = l^/N, where l^ = D C [/2H 2 and D c [ is the 
classical diffusion constant 

D = Uk 2 [1-2J 2 (K)(1-J 2 (K))}, if K>4.5 
cl [0.30(K- K cr f, if Kcr < K < 4.5 ' 

where K cr ~ 0.9716 and J 2 {K) the Bessel function. The localization parameter 
is then defined as (3\ oc — d/N. The scaling law we used is (3\ QC {x) — 7x/(l +72), 
where x = A and 7 4.2 which is slightly different (but in agreement) to 
the one proposed in [7], where x = k 2 /N and K — 5. In Fig. [^fa), we corn- 




Fig. 2. (a) The fit parameter /3br as a function of fi\ oc for 161 x 398 elements, K = 
7(+), 14(x),20(*), 30(D), 35(B) for a wide range of k values, (b) The parameter /3i oc 
vs. A where the scaling law (see text) is shown with the black line. 



pare /?br repulsion parameter of the Pbr(s) with the localization parameter [3\ oc 
through the localization length d =< Ih >■ For the numerical calculations and 
results regarding the spacing distributions P(s) for the quasienergies, we have 
considered a wide range of the quantum perturbation parameter k keeping the 
classical parameter fixed (classically always fully chaotic). In order to amelio- 
rate the statistics, we considered a sample of 161 matrices U nm of size N = 398 
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(« 64,000 elements), in a similar manner as e.g. in [5]) with slightly different 
values of k (Ak = ±0.00125 < k). 

5 Summary 

We studied aspects of dynamical localization in the kicked rotator, following 
[4 5 7], and largely confirm these results. We here considered the case with 
K > 7, where the dynamics is already fully chaotic (ergodic). The fractional 
power law level repulsion is clearly manifested, and globally the level spacing 
distribution is very well described by the Brody or by the Izrailev distribu- 
tion, with a clear and systematic (although not very large) trend towards Brody 
rather than Izrailev. We show that the scaling law (f3\ oc vs. A) exists, but only 
as a first order approximation, as we see some scattering of data around the 
scaling curve. It seems that with increasing dimension of the matrices the scal- 
ing curve asymptotes to the limiting curve with only statistical scattering of 
the data points left. Further research confirms that a similar picture describing 
the dynamical localization applies also in time-independent systems, like e.g. 
billiards [8]. 
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